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^ ! Abstract. 

In the framework of the theory of open systems based on completely positive 
. quantum dynamical semigroups, we give a description of the continuous- variable 

entanglement for a system consisting of two noninteracting modes embedded in a 
thermal environment. By using the Peres-Simon necessary and sufficient criterion for 
, separability of two-mode Gaussian states, we describe the evolution of entanglement 

^ , ' in terms of the covariance matrix for Gaussian input states. For all values of the 

temperature of the thermal reservoir, an initial separable Gaussian state remains 
separable for all times. In the case of an entangled initial Gaussian state, entanglement 
suppression (entanglement sudden death) takes place, for non-zero temperatures of the 
. environment. Only for a zero temperature of the thermal bath the initial entangled 

0^ I state remains entangled for finite times. We also show that, independent of its type - 

separable or entangled, the initial state evolves asymptotically to an equilibrium state 
which is always separable. 
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^ ■ 1. Introduction 

H 



In recent years there is an increasing interest in using continuous variable (CV) systems 
in applications of quantum information processing, communication and computation 
The reahzation of quantum information processing tasks depends on the generation and 
manipulation of nonclassical entangled states of CV systems. A full characterization 
of the nonclassical properties of entangled states of CV systems exists, at present, 
only for the class of Gaussian states. In this special case there exist necessary and 
sufficient criteria of entanglement [21 13] and quantitative entanglement measures [H [5] . 
In quantum information theory of CV systems, Gaussian states, in particular two- 
mode Gaussian states, play a key role since they can be easily created and controlled 
experimentally. 

Implementation of quantum communication and computation encounters the 
difficulty that any realistic quantum system cannot be isolated and it always has to 
interact with its environment. Quantum coherence and entanglement of quantum 
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systems are inevitably influenced during their interaction with the external environment. 
As a result of the irreversible and uncontrollable phenomenon of quantum decoherence, 
the purity and entanglement of quantum states are in most cases degraded. Practically, 
compared with the discrete variable entangled states, the CV entangled states may be 
more efficient because they are less affected by decoherence. 

Due to the unavoidable interaction with the environment, any pure quantum 
state evolves into a mixed state and to describe realistically CV quantum information 
processes it is necessary to take decoherence and dissipation into consideration. 
Decoherence and dynamics of quantum entanglement in CV open systems have been 
intensively studied in the last years P[a[Hlinill0llIIllia[T3l[Tl[T5l|16lll71[Tg 
When two systems are immersed in an environment, then, in addition to and at the same 
time with the quantum decoherence phenomenon, the environment can also generate a 
quantum entanglement of the two systems and therefore an additional mechanism to 
correlate them [H [221 |23l [211 125] . 

In this paper we study, in the framework of the theory of open systems based 
on completely positive quantum dynamical semigroups, the dynamics of the CV 
entanglement of two uncoupled modes (two uncoupled harmonic oscillators) interacting 
with a common thermal environment. The initial state of the subsystem is taken of 
Gaussian form and the evolution under the quantum dynamical semigroup assures 
the preservation in time of the Gaussian form of the state. We have studied 
previously [261 [2Z] the evolution of the entanglement of two identical harmonic 
oscillators interacting with a general environment, characterized by general diffusion and 
dissipation coefficients. We obtained that, depending on the values of these coefficients, 
the state keeps for all times its initial type: separable or entangled. In other cases, 
entanglement generation, entanglement sudden death or a periodic collapse and revival 
of entanglement take place. 

The paper is organized as follows. In Sec. 2 we write the Markovian master equation 
in the Heisenberg representation for two uncoupled harmonic oscillators interacting with 
a general environment and the evolution equation for the covariance matrix. For this 
equation we give its general solution, i.e. we derive the variances and covariances of 
coordinates and momenta corresponding to a generic two-mode Gaussian state. By 
using the Peres-Simon necessary and sufficient condition for separability of two-mode 
Gaussian states [21 [28], we investigate in Sec. 3 the dynamics of entanglement for the 
considered subsystem. We show that for all values of the temperature of the thermal 
reservoir, an initial separable Gaussian state remains separable for all times. In the case 
of an entangled initial Gaussian state, entanglement suppression (entanglement sudden 
death) takes place, for non-zero temperatures of the environment. Only for a zero 
temperature of the thermal bath the initial entangled state remains entangled for all 
finite times, but in the limit of infinite time it evolves asymptotically to an equilibrium 
state which is always separable. We analyze also the time evolution of the logarithmic 
negativity, which characterizes the degree of entanglement of the quantum state. A 
summary is given in Sec. 4. 
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2. Equations of motion for two harmonic oscillators 

We study the dynamics of the subsystem composed of two non-interacting oscillators 
in weak interaction with a thermal environment. In the axiomatic formalism based on 
completely positive quantum dynamical semigroups, the irreversible time evolution of an 
open system is described by the following general quantum Markovian master equation 
for an operator A in the Heisenberg representation (f denotes Hermitian conjugation) 
[291 ED]: 

^ = ^[^'^] + ^UvH^^v^] + [y^'^m)- (1) 

Here, H denotes the Hamiltonian of the open system and the operators Vj,Vj , defined 
on the Hilbert space of H, represent the interaction of the open system with the 
environment. 

We are interested in the set of Gaussian states, therefore we introduce such quantum 
dynamical semigroups that preserve this set during time evolution of the system and in 
this case our model represents a Gaussian noise channel. Consequently H is taken to 
be a polynomial of second degree in the coordinates x, y and momenta Px,Py of the two 
quantum oscillators and Vj,Vj are taken polynomials of first degree in these canonical 
observables. Then in the linear space spanned by the coordinates and momenta there 
exist only four linearly independent operators V^=i, 2,3,4 [31] : 

Vj = a^jPx + ayjPy + b^jX + byjy, (2) 

where axj,ayj,bxj,byj are complex coefficients. The Hamiltonian of the two uncoupled 
non-resonant harmonic oscillators of identical mass m and frequencies Ui and U2 is 

H = ^ipl+pl) + ^iu;W + coy). (3) 

The fact that the evolution is given by a dynamical semigroup implies the positivity 
of the matrix formed by the scalar products of the four vectors a^., b^., a^^, b^^, whose 
entries are the components axj,bxj,ayj,byj, respectively. We take this matrix of the 
following form, where all coefficients Dxx, Dxp^,,--- and A are real quantities (we put 
from now on h = 1): 



(4) 



It follows that the principal minors of this matrix are positive or zero. From the Cauchy- 
Schwarz inequality the following relations hold for the coefficients defined in Eq. (jl]): 

D D - > — D D - D"^ > — 

^xx-^^p^p^ ^xp^ — ^ 5 ^yy^PyPv ^ypy — ^ ' 

DxxDyy - D% > 0, D,^,^D,^,^ - Dl^^ > 0, 

DxxD,^,^ - > 0, DyyD,^,^ - > 0. (5) 
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We introduce the following 4x4 bimodal covariance matrix: 



ait) 



( CTxxit) (^xp^t) (^xy(t) 
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We can transform the problem of solving the master equation for the operators in 
Heisenberg representation into a problem of solving first-order in time, coupled linear 
differential equations for the covariance matrix elements. Namely, from Eq. ([T]) we 
obtain the following system of equations for the quantum correlations of the canonical 
observables (T denotes the transposed matrix) [31] : 
da{t) 



dt 



Ya{t) + a{t)Y' +2D, 



where 
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(7) 



(8) 



(9) 



'yPy ^PyPy / 

The time- dependent solution of Eq. ([7]) is given by [31] 

a{t) = M{t)[a{0) - a{oo)]M^{t) + a{oo), (10) 

where the matrix M{t) = exp(yt) has to fulfill the condition lim^^oo M{t) = 0. In order 
that this limit exists, Y must only have eigenvalues with negative real parts. The values 
at infinity are obtained from the equation 

Ya{oo) + a{oo)Y^ = -2D. (11) 



3. Dynamics of two-mode continuous variable entanglement 

A well-known sufficient condition for inseparability is the so-called Peres-Horodecki 
criterion |2S1 which is based on the observation that the non-completely positive 
nature of the partial transposition operation of the density matrix for a bipartite system 
(transposition with respect to degrees of freedom of one subsystem only) may turn 
an inseparable state into a nonphysical state. The signature of this non-physicality, 
and thus of quantum entanglement, is the appearance of a negative eigenvalue in 
the eigenspectrum of the partially transposed density matrix of a bipartite system. 
The characterization of the separability of CV states using second-order moments of 
quadrature operators was given in Refs. [21 [3] . For Gaussian states, whose statistical 
properties are fully characterized by just second-order moments, this criterion was 
proven to be necessary and sufficient: A Gaussian CV state is separable if and only 
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if the partial transpose of its density matrix is non-negative [positive partial transpose 
(PPT) criterion]. 

The two-mode Gaussian state is entirely specified by its covariance matrix ([6]), 
which is a real, symmetric and positive matrix with the following block structure: 

= ( ^ ) . (12) 

where A, B and C are 2x2 Hermitian matrices. A and B denote the symmetric 
covariance matrices for the individual reduced one-mode states, while the matrix C 
contains the cross-correlations between modes. When these correlations have non-zero 
values, then the states with det C > are separable states, but for det C < it may be 
possible that the states are entangled. 

The 4x4 covariance matrix (fT2l) (where all first moments have been set to zero by 
means of local unitary operations which do not affect the entanglement) contains four 
local symplectic invariants in form of the determinants of the block matrices A, B, C and 
covariance matrix a. Based on the above invariants Simon [2] derived a PPT criterion 
for bipartite Gaussian CV states: the necessary and sufficient criterion for separability 
is S(t) > 0, where 

S{t) =detAdet5 + (^- |detC|)^ 

-Tt[AJCJBJC^J] - i (det A + det 5) (13) 

and J is the 2x2 symplectic matrix 

1 
-1 



j-r,i]- (14) 



3.1. Time evolution of entanglement and logarithmic negativity 

We suppose that the asymptotic state of the considered open system is a Gibbs state 
corresponding to two independent quantum harmonic oscillators in thermal equilibrium 
at temperature T. Then the quantum diffusion coefficients have the following form [30] : 

muiDxx = — = — coth ■ ^ 



mwi 2 2kT 

Dxpx ~ ^ypy ~ Dxy = Dp^py = Dxpy = Dyp^ = 0. 

The elements of the covariance matrix depend on Y and D and can be calculated 
from Eqs. ( llOp . f llip . Solving for the time evolution of the covariance matrix elements, 
we can obtain the entanglement dynamics through the computation of the Simon 
criterion. 

For Gaussian states, the measures of entanglement of bipartite systems are based 
on some invariants constructed from the elements of the covariance matrix [TJ [TTl |33] . 
In order to quantify the degrees of entanglement of the infinite-dimensional bipartite 
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system states of the two oscillators it is suitable to use the logarithmic negativity. For 
a Gaussian density operator, the logarithmic negativity is completely defined by the 
symplectic spectrum of the partial transpose of the covariance matrix. It is given by 
Ej^ = — log2 2i>_, where z/_ is the smallest of the two symplectic eigenvalues of the 
partial transpose a of the 2-mode covariance matrix a : 



A^V^^-4det(T (16) 

and A is the symplectic invariant (seralian), given by A = det A + det B — 2 det C. 
In our model, the logarithmic negativity is calculated as 
1 



EN{t) 



log2[4/(a(t))]. 



(17) 



where 



f{(T{t)) = -(det A + det 5) 



detC 



- (det A + det 5) -detC 



1/2 



det a{t) 



It determines the strength of entanglement for E]^(t) > 0, and if EN{t) < 0, then the 
state is separable. 

In the following, we analyze the dependence of the Simon function S{t) and of the 
logarithmic negativity E]sf{t) on time t and temperature T in the case of a thermal bath, 
when the diffusion coefficients are given by Eqs. (fT5|) . We consider two types of the 
initial Gaussian state: 1) separable and 2) entangled. 

1) In Figure 1 we represent the dependence of the function S{t) on time t and 
temperature T for a separable initial Gaussian state, with the two modes initially 
prepared in their single-mode squeezed states (unimodal squeezed state). Therefore 
the initial covariance matrix is taken of the form 

/ coshr sinhr \ 



.(0) = ^ 



(19) 



sinhr coshr 
coshr sinhr 

\ sinh r cosh r / 

where r denotes the squeezing parameter. We notice that S{t) becomes strictly positive 
after the initial moment of time («S'(0) = 0), so that the initial separable state remains 
separable for all values of the temperature T and for all times. 

2) The evolution of an entangled initial state is illustrated in Figures 2 and 3, where 
we represent the dependence of the function S(t) and logarithmic negativity E^it) on 
time t and temperature T for an entangled initial Gaussian state, taken of the form of 
a two-mode vacuum squeezed state, with the initial covariance matrix given by 

/ coshr sinhr \ 



-(0) = i 





sinhr 




coshr 


— sinh r 





coshr 




— sinh r 


cosh r J 



(20) 
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Figure 1. Separability function S versus time t and environment temperature T for 
A = 0.1, wi = 1,^2 = 3 and separable initial uni- modal squeezed state with squeezing 
parameter: a) r = 1/2; b) r = 2. We take ni = h = k = 1. 



We observe that for a non-zero temperature T, at certain finite moment of time, which 
depends on T, S{t) becomes zero and therefore the state becomes separable. This is 
the so-called phenomenon of entanglement sudden death. It corresponds to the finite 
moment of time when the logarithmic negativity becomes zero. This phenomenon is 
in contrast to the quantum decoherence, during which the loss of quantum coherence 
is usually gradual [I9l |3l]. For T = 0, S(t) remains strictly negative for finite times 
and tends asymptotically to for t oo. Therefore, only for zero temperature of the 
thermal bath the initial entangled state remains entangled for finite times and this state 
tends asymptotically to a separable one for infinitely large times. We notice also that the 
dissipation favorizes the phenomenon of entanglement sudden death - with increasing 
the dissipation parameter A, the entanglement sudden death happens earlier. 

The dynamics of entanglement of the two oscillators depends strongly on the initial 
states and the coefficients describing the interaction of the system with the thermal 
environment (dissipation constant and temperature). As expected, the logarithmic 
negativity has a behaviour similar to that one of the Simon function in what concerns 
the characteristics of the state of being separable or entangled [26], ETJ [3ll 135] . 

3.2. Asymptotic entanglement 

On general grounds, one expects that the effects of decoherence is dominant in the long- 
time regime, so that no quantum correlations (entanglement) is expected to be left at 
infinity. Indeed, using the diffusion coefficients given by Eqs. (|T5|) . we obtain from Eq. 
( ITTl) the following elements of the asymptotic matrices A{oo) and B{oo) : 
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Figure 2. Separability function S versus time t and environment temperature T for 
A — 0.1, wi ~ 1,^2 = 3 and entangled initial vacuum squeezed state with squeezing 
parameter: a) r = 1/2; b) r = 2. We take m = h = k = 1. 



muJ2ayy{oo) = ^PyPy^^^ = 1 coth cryp^{oo) = (21) 

and of the entanglement matrix C(oo) : 

a^y{oo) = (T^p^(oo) = (Typ^ioo) = (Tp^p, (oo) = 0. (22) 
Then the Simon expression (fT3l) takes the following form in the limit of large times, 

S(oo) = i(coth^|^-l)(coth^||,-l). (23) 

and, correspondingly, the equilibrium asymptotic state is always separable in the case 
of two non-interacting harmonic oscillators immersed in a common thermal reservoir. 

In Refs. [26l ETJ Ell EHl |36] we described the dependence of the logarithmic 
negativity E^it) on time and mixed diffusion coefficient for two harmonic oscillators 
interacting with a general environment. In the present case of a thermal bath, the 
asymptotic logarithmic negativity is given by (for uoi < 002) 

E^(oo) = -log2Coth^. (24) 

It depends only on the temperature, and does not depend on the initial Gaussian state. 
Eisf {00) < for T 7^ and En{oo) = for T = 0, and this confirms the previous 
statement that the asymptotic state is always separable. 



4. Summary 

In the framework of the theory of open quantum systems based on completely 
positive quantum dynamical semigroups, we investigated the Markovian dynamics of 
the quantum entanglement for a subsystem composed of two noninteracting modes 
embedded in a thermal bath. We have presented and discussed the infiuence of the 
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Figure 3. Logarithmic negativity En versus time t and environment temperature 
T for ui — \,uj2 — 3 and entangled initial vacuum squeezed state with squeezing 
parameter: a) A = 0.1, r = 1/2; b) A = 0.1, r = 2; c) A = 0.3, r = 1/2; d) A 0.3, r = 2. 
We take m — h = k = 1. 

environment on the entanglement dynamics for different initial states. By using the 
Peres-Simon necessary and sufficient condition for separability of two-mode Gaussian 
states, we have described the evolution of entanglement in terms of the covariance matrix 
for Gaussian input states, for the case when the asymptotic state of the considered open 
system is a Gibbs state corresponding to two independent quantum harmonic oscillators 
in thermal equilibrium. The dynamics of the quantum entanglement strongly depends 
on the initial states and the parameters characterizing the environment (dissipation 
coefficient and temperature). For all values of the temperature of the thermal reservoir, 
an initial separable Gaussian state remains separable for all times. In the case of an 
entangled initial Gaussian state, entanglement suppression (entanglement sudden death) 
takes place for non-zero temperatures of the environment. Only for a zero temperature 
of the thermal bath the initial entangled state remains entangled for finite times, but in 
the limit of infinite time it evolves asymptotically to an equilibrium state which is always 
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separable. The time when the entanglement is suppressed, decreases with increasing the 
temperature and dissipation. We described also the time evolution of the logarithmic 
negativity, which characterizes the degree of entanglement of the quantum state, and it 
confirms the fact that the asymptotic equilibrium state is always separable. 
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